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Abstract 

The evolution of star clusters is determined by several internal and external processes. Here 
we focus on two dominant internal effects, namely energy exchange between stars through 
close encounters (two-body relaxation) and mass-loss of the member stars through stellar 
winds and supernovae explosions. Despite the fact that the former operates on the relaxation 
timescale of the cluster and the latter on a stellar evolution timescale, these processes work 
together in driving a nearly self-similar expansion, without forming (hard) binaries. Low- 
mass clusters expand more, such that after some time the radii of clusters depend very little 
on their masses, even if all clusters have the same (surface) density initially. Throughout it 
is assumed that star clusters are in virial equilibrium and well within their tidal boundary 
shortly after formation, motivated by observations of young (~few Myrs) clusters. We start 
with a discussion on how star clusters can be distinguished from (unbound) associations at 
these young ages. 



1 Introduction 



Ambartsumian (1947) introduced the term association in reference to loose agglomerates and 



he pointed out in subsequent studies that it is unlikely that they are bound by their own 



gravity (see also Blaauw 1964). It is often stated that the majority of stars form in star 



clusters and that there is a high rate of early cluster disruption (e.g. Lada &; Lada 2003). 
In this view associations are clusters that have expanded. But if the star formation process 
is hierarchical then only a small fraction of the newborn stars reside in agglomerates that 



satisfy the conditions necessary to be bound by self-gravity at formation (e.g. Kravtsov & 



Gnedin 2005 Elmegreen 2008 Bressert et al. 2010). When observational samples of star 



clusters are used to support either one of the above scenarios it is vital to know how star 



clusters are separated from associations (Bastian et al. 2011). Here we provide a definition 



of the distinction between these two classes of stellar agglomerates (§ (2J). In § [3] we present 
results of A r -body simulations of star clusters including the effect of stellar evolution. 
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2 The distinction between star clusters and associations 



To illustrate our case we use the recent literature compilation of young massive clusters 



and associations of Portegies Zwart et al. (2010, hereafter PZMG10). This sample consists 



of stellar agglomerates for which a value of the half-light radius r e g, mass M, and age t 
are available in literature. The sample contains 105 agglomerates with M > 10 4 M and 
t < 100 Myr in nearby (< lOMpc) galaxies, including the Milky Way. The ratio of the age 
over the crossing time, II = t/t CT , can be used to separate star clusters from associations 
(PZMG10; Gieles & Portegies Zwart 2011). Objects that are older than their crossing time 



(n > 1) are most likely bound star clusters, whereas objects with II < 1 are expected to be 
unbound associations. The crossing time is defined in terms of empirical cluster parameters 
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where G ~ 4.5 x 10 _3 pc 3 Mq 1 Myr -2 is the gravitational constant. Equation |l| is valid for 
systems in virial equilibrium, because the formal definition includes the root-mean square 
velocity dispersion: t CI oc r e s/(T and in virial equilibrium a oc \J M/r c g. However, a values 
are available for fewer agglomerates and at young ages the empirically determined a can be 



higher because of orbital motions of multiples (Gieles et al. 2010b). Surprisingly, the more 



convenient equation ([!]) fascilites in making the distinction between bound and unbound 
objects. Because super-virial associations expand with a (roughly) constant velocity, their 
crossing evolves as t cr oc r e g oc t. Using equation ([!]) for super-virial objects, therefore, 
overestimates t CT and underestimates EE, thereby pushing them more into the unbound regime. 

In Fig. [T] we show the cumulative distribution of II values of all objects in different 
age bins. The top panel shows that the youngest age bin is a continuous distribution from 
associations with II ~ 0.03 to star clusters with II ~ 10. There seems not to be a distinct mode 
of star cluster formation, but rather a smooth transition from star clusters to associations. 



Bressert et al. (2010) come to a similar conclusion based on the surface density distribution 



of young stellar objects in the solar neighbourhood 

The bottom panel shows that the oldest agglomerates all have II > 1 and these are, 
therefore, star clusters. The intermediate age curves contain both associations and star 
clusters. If we interpret the curves for the different age bins as an evolutionary sequence then 
a distinct gap develops between star clusters and associations around ~ 10 Myr at a value of 
II ~ 1. At older ages an observer should be able to make an unambiguous distinction between 
an (unbound) association and a (bound) star cluster using this straight-forward method. 

For the youngest (continuous) distribution a useful first order separation can still be 
made at a value of IT = 1, as can be noted from the labels of several well known star clusters 
and associations. Independent confirmation comes from recent determinations of velocity 
dispersions of (resolved) young (~ 1 — 2 Myr) star clusters. For NGC 3603 (Rochau et al. 



2010), Westerlund 1 (jMengel & Tacconi-Garman 2007, and M. Cottaar in this volume) and 



R136 in 30 Doradus (Bosch et al. 2009, Henault-Brunet in prep.) it was found that the 
dynamical mass estimates agree very well with the photometric masses, suggesting that these 
clusters are in virial equilibrium and stable (i.e. bound). This is also what their II values 
suggest (see Fig. [I]). In the next section we consider the dynamical evolution of such star 
clusters. 
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Figure 1: Cumulative distribution of II values for nearby (< 10 Mpc) agglomerates in different 
age bins. The vertical dashed line indicates the value II = 1 that separates star clusters from 



associations. From Gieles & Portegies Zwart (2011) 



3 iV-body simulations of star clusters 

We want to understand the evolution of a stellar cluster with a realistic stellar mass function 
in which the stars evolve and lose mass in time ( Gieles et al.||2010a hereafter G10). Because 
most young star clusters are very deeply embedded in their tidal limit (see e.g. Harayama 



et al.|2008 and Gennaro et al.|20lT for the cases of NGC 3603 and Westerlund 1, respectively) 
we ignore the effect of a tidal cut-off. The transition from expansion-dominated evolution to 
tidally limited evolution is considered in more detail elsewhere ( Gieles et al.||20lT ). 

We model clusters with initial half-mass relaxation timescales ranging from t r ho — 
1 Myr to t r ho — 4 Gyr by considering different combinations of [N, log ph] at the start of 
the simulation. The stars follow a Kroupa (2001) initial mass function between 0.1 
and 100 Mq. We use the kira integrator and the stellar evolution package SeBa for solar 
metallicity (Portegies Zwart et al. 2001). The retention fraction of black holes and neutron 



stars was set to zero. 

In Fig. [2] we show how r^/r^Q increases with time as a function of £ r ho- The asymptotic 
behaviour of these results can easily be understood by considering the extremes. Clusters that 
are dynamically young (low i/t r ho) expand adiabatically in order to retain virial equilibrium 
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Figure 2: Expansion from the iV-body runs including the effect of stellar evolution. The full 
lines show a functional description that is used to evolve the inital mass-radius relation in 
Figs. [3] and [| ( see G10 for more details). 

after stellar mass loss. This adiabatic expansion is slow in time and gives a maximum increase 
of ^h/Vho — 2 after a Hubble time. At the other extreme we have clusters that are dynamically 
old (high i/i r ho) an d they expand in a way that is comparable to what is found for equal-mass 
clusters in the sense that all clusters evolve towards the same t r h- 

We conclude that the energy that is released to the core of the cluster as a result of 
stellar mass-loss acts as a central energy source that is in 'balance' with the energy that is 
driven outwards by 2-body relaxation, without the need of forming binaries in the core. One 
of the consequences of this interplay between stellar evolution and dynamical relaxation is 
that there is no sharp transition between a stellar evolution dominated phase and a relaxation 
dominated phase. 

4 Comparing the theory to observations 

With the expression for the evolution of the radius as a function of i r ho at hand (equation 6 
in G10) we can easily calculate the evolution of any initial mass-radius relation. We first 
apply the results to the sample of young objects discussed in the previous section. 

In Fig. [3] we show mass-radius diagrams for all objects for two age ranges: < 10 Myr 
(left panel) and 10 — 100 Myr (right panel). The dashed lines in the left panel are lines of 
constant (arbitrary) p^. In the right panel we have evolved these lines to an age of 100 Myr. 
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Figure 3: Mass-radius relation for all clusters (filled symbols) and associations (open symbols) 
from PZMG10, using = (4/3)r e g. Lines of constant half-mass density, = 3M/(87rr^), 
are over-plotted. The clusters are subdivided into two groups: younger than 10 Myr (left) 
and between 10 Myr and 100 Myr (right). There is a lack of compact (< 2pc) and low-mass 
(< 1O 4 M0) clusters in the older bin due to dynamical expansion. 



These lines are lower limits to the radius at the left-side of this diagram. There are indeed 
no compact (< 2pc) low-mass (< 10 Mq) clusters in this age bin. 

Next we apply our result to the mass-radius relation of old stellar systems in the mass 
range ~ 10 4 - 1O 8 M . In F ig, g we show how a Faber-Jackson type initial mass-radius 



relation (Ha§egan et al. 2005) evolves in time together with data points of (old) globular 



clusters and ultra-compact dwarfs (UCDs) in different galaxies that cover the mass regime 
we are interested in. For high i/i r ho the radius is set by Mq, independent of rho, while for low 
t/trho we are seeing roughly the initial mass-radius relation. At an age of lOGyr the break 
between the two regimes occurs at Mq ~ 10 6 M© and at that age systems with this mass 
have i r h/i — 0.8. Mieske et al. (2008) noticed already that the break occurs at systems with 



i r h roughly equal to a Hubble time. Here we give a quantitative explanation for it. 

We conclude that the evolution of almost all (old) globular clusters is 'balanced', in 
the sense that the rate of central energy production equals the flow of energy due to 2-body 
relaxation^ An important property of this 'balanced' evolution is that the half-mass radius 
is independent of its initial value and is a function of the number of stars and the age only. It 
is therefore not possible to infer the initial mass-radius relation of globular clusters, and we 
can only conclude that the present day properties are consistent with the hypothesis that all 
hot stellar systems formed with the same mass-radius relation and that globular clusters have 
moved away from this relation because of a Hubble time of stellar and dynamical evolution. 



We deliberately avoid the term 'post-collapse' evolution because core-collapse is not required for the 
evolution to be 'balanced', see the discussion in Gieles et al. (20111 
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Figure 4: Mass-radius values for hot stellar systems (see G10 for more details on the data). 
The lines show the evolution of the mass-radius relation. The break at ~ 10 6 Mq at 10 Gyr 
is because lower mass objects have expanded. 
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